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Abstract 

In this paper, we study the behavior of perfect fluid and massless 
scalar field for homogeneous and anisotropic Bianchi type I universe 
model in f(R,T) gravity, where R is the Ricci scalar and T is the 
trace of the energy-momentum tensor. We assume the variation law 
of mean Hubble parameter to obtain exact solutions of the modified 
field equations. The physical and kinematical quantities are discussed 
for both models in future evolution of the universe. We check the 
validity of null energy condition and conclude that our perfect fluid 
solution can behave like phantom model. Finally, we find that perfect 
fluid solutions correspond to massless scalar field models. 

Keywords: f(R, T) gravity; Dark energy; Exact solutions; Massless scalar 
field. 
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1 Introduction 

Over the last decade, the most significant progress in astrophysics and cos- 
mology is the observational evidence that the present universe is undergo- 
ing a phase of accelerated expansion. Observations from supernova type la 
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(SNela) 1 ), cosmic microwave background (CMB) anisotropics 2 \ large scale 
structure 3 -* , baryon acoustic oscillations 4 -* and weak lensing 5 ) show that most 
of the cosmic energy density is dominated by exotic energy source known as 
dark energy (DE). The DE is said to be responsible for this cosmic accelera- 
tion but its nature is still an important challenge. There have been proposed 
many candidates to address this issue. The cosmological constant or vacuum 
energy is the simplest candidate which appears to fit the observational data. 
However, despite of its success the ACDM (A-cold dark matter) model faces 
serious fine-tuning and coincidence problems 6 -*. 

Currently there are two different approaches to address the cosmic accel- 
eration issue. One approach is to introduce various scalar fields of matter in 
Einstein gravity such as quintessence, phantom fields, tachyon field, Chaply- 
gin gas 7 -* etc and also cosmic fluids with anisotropic equation of state (EoS) 8 ). 
The other approach is based on modification of the Einstein-Hilbert action 
to get alternative theories of gravity such as f(R) gravity 9 ), f(T) gravity 10 ), 
Gauss-Bonnet gravity 11 ). Harko et al. 12 ) presented a new modification of 
Einstein Lagrangian by introducing an arbitrary function of scalar curva- 
ture R and trace of the energy-momentum tensor T. The dependence of T 
may be introduced by exotic imperfect fluids or quantum effects (conformal 
anomaly) . 

The action of f(R,T) theory of gravity is given by 12 ) 



where G — c — 1, g is the determinant of the metric tensor g^ u and C m 
determines matter contents of the universe. The energy-momentum tensor 
of matter is defined as 13 ) 



In fact, this modified gravity is the generalization of f(R) gravity and is 
based on coupling between matter and geometry. The corresponding field 
equations have been derived in metric formalism for several particular cases 
of f(R,T) gravity 12 ). They have also explored possibility of reconstruct- 
ing the FRW cosmologies by an appropriate choice of f(T) for the model 
f(R, T) = R + 2/(T). This model has been used to investigate perfect fluid 
solutions of spatially homogeneous and anisotropic Bianchi models 14 ). In 
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recent papers 15 ^, the cosmological reconstruction of f(R,T) gravity models 
have been investigated. We have also explored the validity of first and second 
laws of thermodynamics in f(R,T) gravity 16 ^. 

The existence of anomalies in CMB still requires an intense debate which 
stimulated our interest in anisotropic Bianchi models. Many authors 17-19 ^ 
have studied spatially homogeneous and anisotropic Bianchi models in the 
context of modified theories of gravity. Sharif and Shamir 17 ) explored vacuum 
and non- vacuum solutions of Bianchi I and V universe models in f(R) gravity. 
Bianchi models being anisotropic are useful to study isotropic behavior of 
the universe with the passage of time. Sharif and Kausar 18 ) investigated 
the isotropic behavior of Bianchi III model in f(R) gravity. Scalar field can 
play a vital role to explain the cosmic acceleration which has widely been 
studied in f(R) gravity 20 \ As f(R,T) involves coupling between matter 
and geometry, so considering scalar field as a source may provide some new 
insights. 

In this work, we study perfect fluid and massless scalar field solutions of 
locally rotationally symmetric (LRS) Bianchi type I (BI) universe in f(R,T) 
gravity. Our aim is to find exact solutions of the field equations in this theory 
and discuss the future evolution of the universe. The paper is organized as 
follows: In section 2, we formulate the field equations for LRS BI universe 
in f(R, T) gravity. Section 3 provides solutions of the field equations and 
investigates physical behavior of models and kinematical parameters. Section 
4 contains solutions of the field equations for massless scalar field. Finally, 
section 5 is devoted to discussion and conclusion of the results. 

2 f(R, T) Gravity and Bianchi I Universe 

Variation of the action (TjQ) with respect to the metric tensor results in the 
following field equations of f(R, T) gravity 

Rpi>fit(R, T) - l -g^f(R, T) + (g, u D - V,V u )f R (R, T) 
= 8vrT^ - f T (R,T)T^ u - f T (R,T)e^, (3) 

where fn(R,T) and fr(R,T) denote derivatives of f(R,T) with respect to 
R and T respectively; □ = V^V^, V M is the covariant derivative associated 
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with the Levi-Civita connection of the metric tensor and is defined by 



Since the field equations depend on M1/ , so each form of C m would generate 
a specific set of the field equations. The choice of f(R,T)=f(R) results in 
the field equations of f(R) gravity. The trace of equation is given by 

Rf R (R, T) + 3D f R (R, T) - 2f(R, T) = 8nT - f T (R, T)T - f T (R, T)9, (5) 

where 9 = 6^. This equation leads to f(R,T) in terms of its derivatives 
and matter contents of the universe as 

f(R, T) = i [Rf R (R, T) + 3D f R (R, T) - 8nT + f T (R, T)T + f T {R, T)Q] . (6) 
Eliminating f(R,T) from Eqs. (f5|) and (151) , we obtain 

(R^ - ^Rg^)f R (R, T) + (i^D - V,V u )f R (R, T) = 8n(T, u 

- \Tg,u) - f T (R,T){T, v - l -Tg, v ) - f T (R,T)(e, u - ~9^). (7) 

The line element for the spatially homogeneous and anisotropic LRS BI 
spacetime is given by 

ds 2 = dt 2 - A 2 (t)dx 2 - B 2 (t)(dy 2 + dz 2 ), (8) 

where the scale factors A and B are functions of cosmic time t only. For 
Ait) = B{t) = a(t), this reduces to the flat FRW spacetime. The energy- 
momentum tensor of perfect fluid is 

T^p = {p + p)u^u u - pg^ u , 

where is the four velocity of the fluid, p and p denote the energy density 
and pressure, respectively. We assume equation of state (EoS) p = up with 
u) being a constant. If u = —1, we have ACDM model, u > —1 represents 
quintessence and u < — 1 denotes phantom era. The matter lagrangian can 
be assumed as C m = —p. Using Eq. (plj, Mi/ becomes 

9^ = -2T^ -pg^. (9) 
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Consequently, the field equations ([7]) lead to 

(iV - ±R 9fa ,)f R (R,T) + {- A g^U - V,V u )f R (R,T) = 8n(T^ 
- ~Tg^) + f T (R,T)(T^ + pg^)-^(p + p)g^. (10) 



We consider the function f(R,T) of the form 15 ) 

f(R,T) = f(R) + \T, 



where A is a constant, f(R) is an arbitrary function of R and T, the trace of 
energy-momentum tensor is a correction to f{R) theory. This choice involves 
explicit coupling of matter and geometry, so it can produce significant results. 
The field equations ffTUl) for BI universe model take the form 



A B 
A +2 B 



AB B 2 
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AB B 2 



F + -F 
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^-lj F -2 F -2ll- 2 B) F ^2^ + A ^ + ^^ 



where F(R) denotes derivative of f(R) with respect to the Ricci scalar R 
and 

U = -2 | - + 2^ + 2~ + IL I • (15) 



I + 2 B +2 lB 



S 2 



Now we define some physical quantities for BI model which are important 
in cosmological observations. The average scale factor, volume, expansion 
and shear scalars are defined as 



V = a s 



AB 2 



i 

2°a fc(T 



ab 



. -| 2 
A _ B 

A~B 



• (16) 



5 



The anisotropy parameter of the expansion is characterized by the mean and 
directional Hubble parameters defined as 

A^(^)\ (IT) 

where 

rr „ < 1 (A B 
H = lna = - - + 2— 
v ' 3 yA B 

is the mean Hubble parameter and Hi(i = 1,2,3) represent the directional 
Hubble parameters on x, y and z axes respectively, and are given by 

H. = 4, H, = H. = ?-. 
x A B 

The anisotropy of the expansion results in isotropic expansion of the universe 
in the limit of A — > 0. The deceleration parameter is defined as 

which can be used to explain the transition from past deceleration to the 
present accelerating epoch 21 ). 



3 Solution of the Field Equations 

To solve the field equation in f(R, T) gravity, we assume the variation law 
of mean Hubble parameter defined by the relation 

H = la~ m = l{AB 2 )- m/3 , I > 0, m ^ 0. (19) 

Berman 22 ) proposed this law for spatially homogeneous and isotropic FRW 
spacetime which yields constant value of the deceleration parameter. In re- 
cent papers 16-19 ^, a similar law is used to generate exact solutions for the 
homogeneous and anisotropic Bianchi models in modified theories of grav- 
ity. Using H and V for the BI model in Eq. (1191) . we obtain two different 
volumetric expansion laws 

V = Cl e m , m = 0, (20) 

V = (mlt + c 2 ) 3/m , m^O, (21) 
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where c\ and c 2 are positive constants. Equation (120]) corresponds to de 
Sitter expansion with the scale factor being increasing function of cosmic 
time as a(t) = a^e 111 , H = I = constant. The second volumetric expansion 
law represents power law model with scale factor a{t) = aot n . If < n < 1, 
then power law solution is accelerating and for n > 1, it exhibits decelerating 
behavior. Subtracting Eq.f lT4"|) from (fT3"|) with some manipulation, it follows 
that 

k 

H x -H y = —, (22) 
where A; is a positive constant. Using Eqs. ffT6]) and f )22|) in ffTT|) . it turns out 

* V. (23) 



■0) \V3VF 
In the following, we discuss above two cases separately. 



3.1 Perfect Fluid Model When m = 

Here, the spatial volume of the universe for exponential expansion is given 
by Eq.(l20l). Using this value of V in Eq.(|22l). we can write the scale factors 
as 

A = cl /3 c 2 3 /3 e lt+ f B = cl /3 cz 1/3 e lt -y ^ dt . (24) 

To find the explicit solution of the field equations, we assume a relation 
between F and a as F ex a nl7 ~ 18 \ which implies that 

F = ae nl \ 

where a is the proportionality constant and n is any arbitrary constant. As 
we are interested to discuss the exponential and power law expansions, so 
it would be useful to assume unknown F in terms of these expansion laws. 
This assists to reconstruct the f(R, T) gravity depending on the choice of 
the scale factor. Using this value of F in Eq. (|24p . we obtain 

, 1/3 2/3 ^- 2* e -("+3)" 1/3 -1/3 lt+ , ,, k ,„ e -("+3)« , . 

For n > —3, we observe that the scale factors A(t) and B(t) have constant 
values at initial epoch which imply that the model has no initial singularity, 
while these diverge in future evolution of the universe. When n < —3, the 
scale factors increase with time and approach to very large values as t — > oo. 
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Figure 1: Plot of A versus cosmic time for different values of n. We set 
I = 0.1, k = 3, and a = 0.05. (Colour online) 



For n = —3, the model represents similar behavior in every direction. The 
directional, mean Hubble parameters and anisotropy parameter of expansion 
turn out to be 



H 
A 



2k 



x -1 + 

6a 



-(n+3)Zt 



Hy = H Z 



k 
( 

3a 



-(n+3)lt 



H = L 



2k 2 
9l 2 a 2 ' 



,-2(n+3)«t 



The mean Hubble parameter is found to be constant whereas the directional 
Hubble parameters are dynamical. For n > —3, H x and H y become constant 
at t = as well as for t — > oo. These parameters vary from H by some 
constant at t = and coincide for later times of the universe. As the constant 
being positive (negative), it would increase (decrease) expansion on the x- 
axis and it decreases (increases) expansion on y and z axes. For n = —3, H x 
will increase from if by a constant factor while parameters H y , H z will 
decrease by a factor The anisotropy parameter of expansion is found to 
be finite for earlier times of the universe and vanishes as t — > oo for n > —3. 
The plot of A versus cosmic time t is shown in Figure 1. 

The deceleration parameter, expansion and shear scalars are given by 



q=-\, 6 = 3l = 3H, 



a 



3a 1 



-2(n+3)lt 



(26) 
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Figure 2: Evolution of NEC for n = 2. The left graph shows that NEC is 
satisfied for a < and it is violated for a > at the right side. We set 
/ = A = 0.1 and k = 3. (Colour online) 

The volume V of the universe is an exponential function which expands 
with the increase in time and becomes infinitely large for later times of the 
universe. Also, the expansion scalar is constant for all times and hence 
the model would favor the uniform expansion. The deceleration parameter 
(q = —1) allows the existence of accelerating model for this case which is in 
agreement with the current observations of SNela and CMB 1-2 ). 

Using Eqs. (I25j) in Eqs. (I12 |) - f ll4p . we obtain the following relation of energy 
density and pressure 

p + p = ~ 1 - \2k 2 e nlt - 2{n+3)lt + 3b ie nlt ] , (27) 
3a(87r + A) L J 

where b\ = n(n — l)l 2 a 2 . This shows that the null energy condition (NEC) 
is violated, i.e., p + p < which implies that ui < — 1. Matter component 
with u < —1 is named as "phantom energy" and is a possible candidate 
of the present accelerated expansion. The phantom regime favors recent 
observational cosmology of accelerated cosmic expansion. The behavior of 
NEC for different values of a is displayed in Figure 2, which shows that NEC 
is violated for positive values of a. Thus, we assume a > for phantom 
universe. Equation ([27)1 implies the following dynamical variables of the 
perfect fluid 
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Figure 3: The left graph shows the behavior of p for — 6 < n < and 
n < —6, while the right graph presents the evolution of p for n > 0. We set 
I — A — 0.1, k = 3 and a = 0.05. (Colour online) 

- j [2k 2 e nlt - 2{n+3)lt + 3b ie nlt ] , (28) 



r 3a(l + w)(87r + A) 

^ 3a(l + w)(87r + A) L 

For the phantom evolution of the universe, p decreases with cosmic time 
and approaches to zero as t — > oo in the range of —6 < n < 0. When n < —6, 
p increases as time goes from zero to infinity and hence diverges. Figure 3 
shows that p decreases for n = —1, —4 and becomes uniform for n = —6. 
However, the value of n = —7 shows increasing p for the future evolution of 
the universe. If n > 0, p decreases with time but for large values of n, it 
shows bouncing behavior as shown in Figure 3. For this model, the scalar 
curvature R and f(R, T) are given by 

R = _ ^ r 18zV + A .2 e -2(«+3)ttl 

3cr L J 

f(R,T) = ™(R + 3l(nH + 3))e"« + ^ 7 ^ + ^ 
71 ' ; 2 V v ;; 6a(87r + A)(l + w) 

x (2A; 2 e- (n+6)H + 36 ie n "). 
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3.2 Perfect Fluid Model When m^O 



For m 7^ 0, the spatial volume is given by Eq. (j21|) and the corresponding 
deceleration parameter is q — m — 1. If ^ > 0, the model represents deceler- 
ating universe whereas q < indicates inflation. To obtain the accelerated 
expansion model, we take m < 1. Solving the field equations ( JT2|) -(JT4 ]) . the 
scale factors are found to be 



n + 3 



A = cj\mlt + C2 ) 1/m e^(™-"-v (mlt+C2)1 ~ 

, n-|-3 

£? = c ~ 1/3 (mlt + ^l^e ^ra ^i " m . (30) 

We discuss the evolution of the scale factors for two cases, i.e., m > n + 3 
and m < n + 3 along with < m < 1. If m > n + 3, the scale factor ^4 
increases with time whereas B tends to zero. For m < n + 3, the behavior 
of scale factors is almost identical provided that n is always greater than —3 
to keep m positive. Substituting Eq. (l30|) in Eq. lflTl) . we get 

2k +3 
H = l{mlt + c 2 )-\ H x = l{mlt + C2)' 1 + —{mlt + c 2 y^ , (31) 

6a 



k 



n + 3 



H y = H z = l{mlt + ca)- 1 - -^Ajnlt + c 2 )~— , (32) 

2fc 2 . , . o (n+3) 

A = 9l^ {mU + C2) ■ (33) 

The Hubble parameters H, H x , H y and H z become constant at the initial 
epoch. As t — > 00, the values of these parameters tend to zero for n > — 3 
and become infinite for n < —3. 

If n < —3, A increases with cosmic time whereas for n > —3, its value 
decreases and may result to isotropic expansion in future evolution of the 
universe (see Figure 4). The expansion and shear scalars turn out to be 



k 2 , ,. . ,_9(=±3) 

3cr 



9 = 3l(mlt + c 2 ) , <T = — imlt + c 2 )~ 2 ~. (34) 



If we replace Eq.flHU]) in Eqs.fllSD-CED, we obtain 



P + P 



-1 



3a(87r + A) 



(n + 6) (n-2m) 



2k 2 imlt + c 2 ) — — + 3b 2 imlt + c 2 )^^~ ,(35) 
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Figure 4: Plot of A versus cosmic time t for different values of n. We set 
I = 0.1, k = c 2 = 3, m = 0.9 and a = 0.05. (Colour online) 



where b 2 = (n(n — 1) — m{n + 2)) la. p and p are obtained as follows 



•1 



P 



3a (1 + u)(8tc + A) 
— w 

3a(l + w)(87r + A) 



2k 2 (mlt + c 2 ) 
2k 2 (mlt + c 2 ) 



(n+6] 



+ 3b 2 (mlt + c 2 



(n+6) 



+ 3b 2 (mlt + c 2 



(n— 2m) 



(n — 2m) 



(36) 
(37) 



Equation (135 j) shows that NEC is violated for the power law expansion 
model. The behavior of NEC is shown in Figures 5-6 which depends on the 
choice of a as well as n. For each value of n except 2 < n < 0, NEC can be 
satisfied for a < but the choice a > does not support it. If 2 < n < 0, the 
constraints to satisfy and violate NEC are interchanged (see Figure 6). For 
m < 1, energy density decreases in the range of — 6 < n < and increases 
with cosmic time t for n > —7. For n > 0, the behavior of p is shown in 
Figure 7. From Eqs.(|6]) and (115]) . the Ricci scalar i? and f(R,T) are given 
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Figure 5: Behavior of NEC versus a for n = 3. The left part shows that 
NEC is satisfied for a < 0, while it is violated for a > shown on the right 
side. We set I = X = 0.1, k = c 2 = 3 and m = 0.9. (Colour online) 




Figure 6: This figure is plotted for n = 2. The left part shows that NEC is 
violated for a < 0, whereas NEC is satisfied for a > shown on right side. 
We set / = A = 0.1, k = c 2 = 3 and m = 0.9. (Colour online) 
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Figure 7: Plot of p versus cosmic time t for n > 0. We set / 
k = c 2 = 3, m = 0.9 and a = 0.05. (Colour online) 



A = 0.1, 



by 



R 



3a 2 

Oi 



9l 2 a 2 (m - 2)(mlt + c 2 y 2 - k 2 (mlt + c 2 )~ 2( ™ } 

2 



/(i?,T) = - [R(mlt + c 2 )™ + 3n(n-m)l 2 (mlt + c 2 )™ 

, x™ n 87r(l — 3a;) + A(l — u;) 9 , , . n+e 
+ QHmlt + c?)" 1 ] + ; /o _ , — - (2k 2 (mlt + c 2 ) 



6a(8yr + A)(l + w) 



+ 36 2 (m/t + c 2 )-" 2 ). 



4 Massless Scalar Field Models 

The Lagrangian for massless scalar field 4> is given by 20 ) 
and the corresponding energy-momentum tensor is 



(38) 



= d^(f)d v (f) - -g^dyipd 1 ^. 



(39) 



14 




Figure 8: Evolution of versus cosmic time t for m = and different values 
of n: solid (black) n — 1; dashed (red), n = 0; dahsed(blue), n — — 1. We set 
/ = A = 0.1, ft = 3 and a = 0.05. (Colour online) 



Here, T^ v represents stiff matter with EoS uj^ = 1. Using Eqs.fJZJ), (JHJ) and 
(|39|) . we obtain the following field equations for massless scalar field 

4 - 2* + 2^ - 2* V - If + ( 3 4 - |V = 1(8, + 
t4 B AB B 2 J 2 ^2A i3 y 2 V ^ ' v ; 




4.1 Scalar Field Models When m = and m/0 

The field equations with massless scalar field are similar to the perfect fluid 
case with Ptf, + p<f> = (ft 2 . Hence, we obtain the same results for the scale 
factors and other physical parameters. Substituting Eq.()25l) in (H0l -( H2l) . we 
obtain the time derivative of scalar field 6 as 



± / 1 \ 2k 2 e nlt-2{n+Z)lt + 3fe nit] (43) 

w 3a(8vT + A) [ J K ' 
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Figure 9: Evolution of <fi versus cosmic time t for m^fl and different values 
of n: solid (black) n — 1; dashed (red), n = 0; dahsed(blue), n = — 1. We set 
Z = A = 0.1, k = c 2 = 3, m = 0.9 and a = 0.05. (Colour online) 



Using Eqs.@25J and ([43]) in (J6J), it follows that 

47T 



/(*.n = f(*+ 3))**-. 



(2A; 2 e- (n+6)/ * + 3b ie n "). (44) 



The behavior of for exponential expansion is shown in Figure 8. If we 
solve Eqs. (j40p - ()42p for the case m / 0, we get the similar solutions as given 
in section 2.2. The expression of <ft is obtained as follows 



1 



3a(8vr + A) 



2k 2 (mlt + c 2 ) ~ + 3b 2 (mlt + c 2 ) ™ 



(45) 



Evolution of versus cosmic time t for different values of n is shown in Figure 
9. Substituting Eqs.JSU and (BHD in E q-©, we have 



f(R, T) = | [i2(mft + c 2 ) « + 3n{n - m)l 2 {mlt + c 2 )™~ 2 + 91 (46) 



x (mlt + c 2 )™ x ] 



47T 



3a(87r + A) 



n+G 



2k\mlt + c 2 ) ~ + 3b 2 (mlt + c 2 ) 
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5 Discussion and Conclusion 



The issue of accelerated expansion of the universe can be explained by taking 
into account the modified theories of gravity such as f(R,T) gravity 23 ^. In 
f(R,T) gravity, cosmic acceleration may result not only due to geometric 
contribution to the matter but it also depends on matter contents of the 
universe. The coupling between matter and geometry in this gravity results 
in nongeodesic motion of test particles and an extra acceleration is always 
present. This theory can be applied to explore several issues of current inter- 
est and may lead to some significant results as compared to other modified 
theories. 

The homogeneous and anisotropic Bianchi universe models with perfect 
fluid have been investigated in the context of f(R,T) gravity 14 ). The exact 
solution of the field equations are obtained for the particular choice f(R, T) = 
R + 2/(T) with f(T) = XT. The results of this formulation are very similar 
to that in general relativity and cannot imply the reconstruction of f(R,T) 
gravity 14 -*. Houndjo 15 ) used a more general form f(R,T) = f(R) + AT to 
reconstruct the f(R, T) gravity from holographic DE numerically. In present 
work, we have employed this choice to reconstruct some explicit models of 
f(R,T) gravity for BI universe. We have presented the phantom evolution 
of the universe by examining the NEC. 

The exact solutions of the modified field equations are obtained for the 
spatially homogeneous and anisotropic LRS BI universe with perfect fluid 
and massless scalar field. The law of variation of mean Hubble parameter is 
assumed which implies two cosmological models for m = and m^O. These 
models represent the accelerated expansion of the universe which supports 
the observations of WMAP data and SNela 1-2 ). We have presented physical 
properties of the models as well as kinematical parameters. In the following, 
we summarize the results for these two models. 

• Model for V = Cie m 

For exponential expansion model, the accelerated expansion of the uni- 
verse may occur since q — — 1. The kinematical paramters have been dis- 
cussed for two cases n > — 3 and n < —3. The expansion scalar is constant, 
while the Ricci scalar approaches to constant value as t — > oo for n > — 3 
and take infinitely large values for n < —3. The anisotropy parameter of 
expansion depends on time and vanishes in future evolution for n > —3. If 
a > 0, NEC is violated. In our discussion, we have considered a particular 
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model of f(R, T) gravity, f(R, T) = f(R) + AT. We are not able to find the 
explicit function of f(R,T) by using Eq.flB]). For A = 0, we develop f(R) in 
terms of R and hence the function f(R, T) as 



f(R) 



(R + 12H 2 ) + ^{3a 2 H((n 2 - A)H + 3) + V"^M 



e nlt 



l + ' 2a L ' 1 + oj 

which can be expressed as 

f(Rx) = constt x R x mi + const 2 x R x m \ 

where R x — R+ 12H 2 , mi = j^+g) an d m 2 — m i ~ 1- The models of f(R\) 
depending on n are shown in table 1. We find that: 

Table 1: Models of f(Ri) corresponding to n. 



n 


f(Ri) 


n = 


Rx + Rl, R^ = const 


n = -2 


R\ + R\ 


n = —6 




n = — 4 


R x + R x ~ 


n = ~f 


R'i + R( 


n = ~f 


+ i?r a 


n = -§ 


i?f + R? 


n = -| 


Rj + 



For n = 0, /(i?) represents the ACDM model, i.e., = i? + A. If we 

put constant = 0, then /(/?, T) is of the form f(R, T) = R + T. The most 
famous Starobinsky's model 24 ), f(R) = R + aR 2 is achieved for n = —2 and 
the corresponding f(R, T) function is f(R, T) = R + aR 2 + T. For n = —6, 
f(R, T) can be presented as f(R, T) — j^+T. 

For massless scalar field (m = 0), we have found similar results for scale 
factors as in perfect fluid. The expression of f(R) is 



f(R) 



a{R + 12H 2 ) + —{3a 2 H((n 2 - A)H + 3) - bj} 
2a 



i.e. 



f{R\) = const 3 x i?! mi + const A . x R{ 
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• Model for V = (mlt + c 2 ) 3/m 

For m/0, the deceleration parameter is q — m — 1, which leads to the 
accelerating universe model for < m < 1 and if m > l(q > 0), the model 
represents decelerating phase of the universe. The evolution of the scale 
factors is discussed for two cases m > n + 3 and m < n + 3 with < m < 
1. The anisotropy parameter of expansion increases for n < — 3, whereas 
it may result in isotropic expansion in future evolution of the universe for 
n > — 3. The Hubble parameter, expansion scalar and shear scalar approach 
to constant at earlier times of the universe and approach to zero as t — > oo. 
The scalar curvature R becomes constant as t — > oo for n > —3, whereas it 
diverges for n < — 3. For A = 0, we have 



f(R) = (mlt + 02) 



n/m 



J^L(R _ 6 (m - 2)H 2 ) + -H{3H(n 2 



X 3c^ 

+ m ~ 4 )+ 2/2 n 1 > # + 9} 
crr(l + 

which results in 

/(i? 2 ) = const 5 x .Rs™ 1 + const 6 x R 2 m \ 
where R 2 = R — 6(m — 2)H 2 . In case of massless scalar field, f(R) is 

f(R) = (mlt + c 2 ) n/m \a(R-6(m-2)H 2 ) + ^H{3H(n 2 + m 
- 4) + ^|| + 9} 

i.e., f(R 2 ) = const 7 x R 2 mi + constg x i? 2 ™ 2 - The scalar field is found to 
be decreasing function of cosmic time (see Figures 8, 9) 19 ^. 

We have seen that all f(R) represent identical behavior with different 
constraints. The NEC is found to be violated for both models m = and 
m^O which results in phantom evolution. For uj < — 1, energy density is 
found to be positive and pressure is negative. Thus, our solutions for perfect 
fluid represent the phantom era of DE which are consistent with present 
observations of WMAP5 25 ). The isotropic behavior of models is observed 
for future evolution. It can be concluded that solutions of massless scalar 
field in f(R, T) gravity can be recovered for both exponential and power law 
expansion models from phantom solutions of perfect fluid if cu — 1. 
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